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Abstract 

Several examples are known where quantum gravity effects resolve the classical 
big bang singularity by a bounce. The most detailed analysis has probably occurred 
for loop quantum cosmology of isotropic models sourced by a free, massless scalar. 
Once a bounce has been realized under fairly general conditions, the central questions 



OJT)| are how strongly quantum it behaves, what influence quantum effects can have on 

its appearance, and what quantum space-time beyond the bounce may look like. 
This, then, has to be taken into account for effective equations which describe the 
evolution properly and can be used for further phenomenological investigations. Here, 
we provide the first analysis with interacting matter with new effective equations valid 
for weak self-interactions or small masses. They differ from the free scalar equations 
by crucial terms and have an important influence on the bounce and the space-time 
around it. Especially the role of squeezed states, which have often been overlooked 
in this context, is highlighted. The presence of a bounce is proven for uncorrelated 
states, but as squeezing is a dynamical property and may change in time, further 
work is required for a general conclusion. 

X 

1 Introduction 

Classical cosmology as a description of the early universe is plagued by the big bang 
singularity. In this extreme high energy and high curvature regime quantum space-time is 
expected to take over the reigns from general relativity in a way which will hopefully cure 
the singularity problem. A common expectation is that the universe, instead of collapsing 
completely, will bounce at a minimum non-zero volume and thus connect to a well-defined 
pre-big bang state. Mediated by quantum gravity, such a picture would exploit the changes 
in forces and the dynamics which quantum physics should imply. If this can be realized, 
a central question remains: Do the quantum forces that have been summoned to fight 
the classical singularity leave space-time otherwise intact, or will space-time remain in a 
fluctuating turmoil long after the immediate danger of divergences has passed? Generally, 
strong quantum effects as they may be necessary to prevent the singularity lead to large 
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fluctuations. Then, a genuine quantum space-time has to be dealt with, which may not be 
near a smooth bouncing effective picture. 

Detailed examples have been provided by loop quantum cosmology [lj, where homo- 
geneous models are governed by a non-singular difference equation for quantum states 
[21 El IH [5] . The general difference equation, however, does not provide direct insights into 
intuitive geometric properties around the classical singularity; in general the transition may 
not be semiclassical at all. But there is a class of models which are more easily accessible: 
isotropic models sourced solely by a free, massless scalar. Initial numerical calculations 
P, [7] have suggested that the singularity is avoided by a surprisingly smooth bounce of a 
wave packet which remains sharply peaked and fluctuated before the bounce very nearly 
as it does afterwards. "Before" and "after" here refers to a global internal time which in 
these models is provided by the value of the scalar. 

Upon closer examination, however, two new key features materialized: 

• In a suitable factor ordering of the Hamiltonian constraint operator, a flat isotropic 
model with a free, massless scalar is a solvable system in which quantum back- 
reaction does not occur [Sj. This is the case even in the bounce phase of the loop 
quantization. Thus, quantum variables such as fluctuations or higher moments of 
a state do, although they are themselves dynamical and change in time, not influ- 
ence the evolution of expectation values. It is this feature which makes the models 
highly controlled, but also very special. In quantum mechanics, for instance, this be- 
havior is realized only for the harmonic oscillator. When interactions, or deviations 
from isotropy, are included quantum back-reaction does result and the evolution will 
become more quantum in the sense that it depends on how the quantum variables 
behave. Quantum correction terms due to coupling between expectation values and 
fluctuations would at first not be large for evolution starting with a semiclassical 
state, as it should be realized at large volume of a universe. However, states spread 
and deform from an initial Gaussian distribution. Accordingly, quantum variables 
generically grow in time and with them quantum corrections. This may take time, 
but can be significant in long-term evolution. Cosmology is a study of the system 
with the longest possible evolution times, and especially regarding the big bang sin- 
gularity the correct question to ask is: What is the behavior of a universe at small 
volume, described by a generic quantum state as it has evolved from a semiclassi- 
cal state at large volume? This question cannot be answered if the behavior of the 
quantum system is known only for semiclassical states and rather short evolution 
times. For all we know, the universe near the big bang can be in a highly quantum 
state. Thus, quantum back-reaction is a serious player in this regime, and one has 
to understand how it may affect the bounce of the free solvable model. 

• Even within the solvable model, fluctuations before and after the bounce can be very 
different from each other [9j [TUJ This asymmetry is not a consequence of quantum 

1 Some of the recent literature on this topic has witnessed a certain amount of reluctance to accept this 
property. For instance, it has been claimed in [llj that fluctuations of a state which is semiclassical at 
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back-reaction but of squeezing. As the precise state is unknown, there is no strong 
control on how it may be squeezed even when it is semiclassical. Accordingly, no 
strong control on fluctuations during and before the big bang exists. Still, for any 
such state there is a bounce in this free model provided that the state is semiclassical 
at one time. 

The formulation in [8] as a solvable model allows the derivation of precise effective 
equations in this case, which easily shed light on the evolution of a quantum space-time and 
its properties. If the model is perturbed by including a small matter potential, perturbative 
effective equations can still be derived in an expansion [12]. In this paper, we use these 
equations to show that the bounce is in fact affected by interactions and determine how this 
is related to the possible squeezing of the state. Our analysis is thus based on a combination 
of the two key results which are realized in isotropic, free scalar models. We provide a new 
effective equation which, unlike previous examples, is valid for a massive or self-interacting 
scalar. It contains additional correction terms which are especially important near a would- 
be bounce. 

2 Loop quantum cosmology 

In loop quantum cosmology, the basic pair of canonical variables is formed by the square 
p = a 2 of the scale factor a and extrinsic curvature, i.e. the time derivative c = 7a. Their 
Poisson bracket is {c, p} = 8^7(^/3. These are isotropic reductions of a densitized triad 
Ef = pbl and the Ashtekar connection A % a = c5 l a and enter the basic operators of a loop 
quantization. Although p can take both signs, signifying the two possible orientations of 
the triad, we will mainly restrict attention to positive values. 

While p, just as the densitized triad of the full theory, can readily be quantized through 
flux operators, there is no operator for the connection components. Instead, holonomies of 
the connection are represented in a well-defined way, which are obtained after integrating 
the connection along curves and exponentiating the result. In an isotropic model, it is 
not obvious what the length of the curve used in the integration should be, thus giving 
rise to quantization ambiguities which can only be resolved in a detailed reduction of the 
isotropic model from the full quantum theory. (See [EH EH ESI ESI E3 ESI EH] for work 
in this direction, which remains incomplete.) What lengtlH is used should depend on the 
underlying state and its scale of discreteness, which in general should be time dependent. 

large volume have to be very nearly symmetric around the bounce. What was actually derived there was 
an upper bound for the difference before and after the bounce of fluctuations of volume relative to the 
total volume. However, the discussion ignored the size of these relative fluctuations and overlooked the fact 
that each of the terms in the difference, not just the difference, is smaller than the upper bound provided. 
Thus, the inequality does not at all restrict the asymmetry. One can easily see that the numerical example 
provided in [11] allows for ratios of relative fluctuations before and after the bounce as large as 10 28 . 
Judging from this estimate alone, a state before the big bang may be highly non-semiclassical even if it 
becomes very semiclassical after the big bang. The detailed bounds provided in [lOj are much sharper, but 
still leave room for asymmetric fluctuations. 

2 The notion of length here only refers to coordinates and so it may seem that it is not covariant. 
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Thus, also the length factor /x of a holonomy exp(z/ic) in an isotropic model could be a 
function of time, or rather of the spatial geometry since this determines the evolving scales. 

A quantum theory with a time dependent number of degrees of freedom is difficult 
to formulate and analyze [20, ED [22], [23] EI] , but key properties can be included in the 
description of an isotropic model. Instead of using a fixed fi G R, we formulate the 
dynamics with holonomies of the form exp(if(p)c) where f(p) takes into account a possible 
dependence of the discreteness scale on the total size. When quantized, the exponential 
of the momentum c becomes a shift operator in p, rather than a derivative operator. 
For f\p) = const, the basic triad variable p would be equidistantly spaced when acting 
with powers of holonomies, but not for a p-dependence of f(p). Of special interest are 
functions f(p) = fop x of power-law form, which can arise from loop quantum gravity for 
— 1/2 < x < [T7j. Several independent phenomeno logical and stability arguments prefer 
a value near the lower bound x ~ — 1/2 [25] 126] |2T] 128] 129] 130] . For arbitrary x, a canonical 
pair is formed by the new variables {f(p)c, V) with 

3» 3a 2 ( 1-x ') 

v '■= o — Tvi — = s — ru~\ W such that {/oP*c, V} = 1 . (1) 

8tt7G(1 - x)f{p) 8tt7G(1 - x)f 

(We have V proportional to the spatial volume for the power- law case f(p) oc p~ x l 2 . In 
this case, the number of vertices in a graph would grow in a way proportional to spatial 
volume: Af(a) = f(a 2 )~ 3 oc a 3 . ) This shows that for a given x, p/f(p) is equidistantly 
spaced upon action of exp(z/(p)c). 

In loop quantum cosmology, we have to quantize the classical Friedmann equation 

2 

a\ 8tcG H mattPr (a) 8nG 

- = -z~P (2) 



a J 3 a 3 

where -f/matter is the Hamiltonian of matter and p its energy density. This is done via the 
Hamiltonian constraint, which is obtained by multiplying the Friedmann equation with a 3 . 
For a free, massless scalar, we have the matter Hamiltonian -ff ma tter = | a_3 P0 with the 
scalar momentum p^. The gravitational part requires the use of c, which in loop quantum 



However, no such problems or even issues with potential gauge artifacts arise as one can see upon closer 
examination [17]. To see this, we have to refine our definition of the basic variables as p = V^ 3 a 2 

1/3 

and c — V "fa where Vq is the coordinate volume of a region chosen for the spatial integration of the 

— 1/3 

symplectic structure and the Hamiltonian |4J. Holonomies then take the form exp(i£oV c) where Iq is 
the coordinate length of an edge as referred to above. This immediately shows that the expression is not 
coordinate dependent. It also does not depend on the choice of Vq since this dependence cancels between 
c and V Q 1 ^ 3 . Instead of referring to a coordinate quantity, the length factor can be interpreted as A/" -1 / 3 
where J\f = Vq/£q is roughly the number of vertices of an underlying spin network state in the volume 
Vo- This immediately shows the dynamical aspect of a pre-factor /(p) in exp(if(p)c), which includes the 
factor Vo, if the graph is being refined during evolution. In this way, Af(p) will be the number of vertices 
of the state at volume p 3 ^ 2 . In this context, notice that a solution to the Hamiltonian constraint equation 
in general is a superposition of many states associated with different graphs. The function Af(V) refers 
to this relationally as the number of vertices of graphs in a decomposition as eigenstates of the volume 
operator. This indicates how f(p) refers to internal time evolution in a constrained theory. 



4 



cosmology can only be formulated by means of holonomies which are periodic functions 
of their arguments. There is no exact correspondence, but all we need is a well-defined 
operator which has ([2j) as its low-curvature limit. Many examples for such operators exist, 
the simplest choice being of the form 

g^(/(p)- J i5(/(p)<:)v5)#^) = Afv^^M) ■ (3) 

(See e.g. [3T] for a discussion of some ambiguities in this case.) We have not yet decided 
on the factor ordering between p and c on the left hand side, and also the right hand 
side is not fully defined as written because p~ 3 / 2 requires an inverse, which does not exist 
for the operator p with its discrete spectrum containing zero in loop quantum cosmology. 
Its quantization thus has to be more indirect, which can be done [32] following general 
constructions of the full theory [33] . (Also here, ambiguities arise [MJ EH] without changing 
the qualitative picture.) 

2.1 Difference equation 

Basic operators of the quantization are the flux operator p\p) = ^7i^£ F p\p) with the Planck 
length lp = y/Gh, whose orthonormal eigenstates \p) with p G M we are using, and the 
holonomy operator exp(ip'c/2)\p) = \p + p'). By expressing the Friedmann equation in 
terms of holonomies, which then become finite shift operators, one thus obtains a difference 
equation for wave functions [251 E] • 

There are several ingredients for an explicit expression of the difference equation: First, 
using holonomies for a given refinement scheme we refer to shift operators exp(if(p)c/2) \p) ~ 
\p + F(p)) of a /i-dependent step-size, where F(p) := f^n^ipp). (This action is not pre- 
cisely written yet due to factor ordering which is discussed below.) The factor */p in (EJ) 
could be quantized directly via p, but is usually written in a way closer to what one has 
in the full theory. The reason is that in an inhomogeneous expression one would instead 
have a factor E^Ej/ yj det(E%)\ which, due to the inverse of the determinant, cannot be 
quantized directly. Following the same steps in an isotropic model, one is led to a quanti- 
zation of y/p of the form if(p)~ 1 exp(if(p)c)[exp(—if(p)c), \p\ 3 ^ 2 ] which also happens to be 
diagonalized by the triad eigenstates \p). (Note, however, that the diagonalization prop- 
erty does not hold in the full theory [33, [38].) Finally, the factors f(p)~ 2 and p~ 3 / 2 have 
to be turned into operators. 

An explicit realization faces two main issues: (i) the non-equidistancy of step-sizes and 
(ii) the factor ordering. The ordering between yfq and f(p)~ 2 sin 2 (/(p)c) can be mimicked 
from the full theory [39, 40J, which leads one to order the quantization of yfp to the right 
to obtain the basic Hamiltonian constraint operator H which may in a second step be 
ordered symmetrically as ^(H + H^). However, the factors of f(p) model the refinement 
behavior of an underlying lattice state, and thus do not occur in the full theory where the 
lattice refinement is automatically realized by elementary holonomies. Thus, the ordering 
in isotropic models is not strictly defined, and one cannot always use guidance from the full 
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theory. In applications, one should therefore use only properties which are insensitive to the 
factor ordering. This is also useful for the first issue: the non-equidistancy of the resulting 
difference equations. As shown in [26], one can always transform a difference equation of 
isotropic models to equidistant form, up to changing the factor ordering. Thus, one can 
use preferred orderings in which an equidistant difference equation results which is much 
easier to analyze. (See [JT] for a numerical procedure to analyze non-equidistant difference 
equations directly.) 

Specifically, we obtain a Hamiltonian constraint operator 

{H-H m )\im) = 3/2 ^ (l M +F( / x)| 3 / 2 -l M -F( M )l 3 / 2 )F( M )- 3 (i M +4F( M ))-2l/x)+l M -4F( M ))) 

^ (4) 
by following the steps detailed in [HH]. A physical state \tp) = ^2 M (0)]/i) has to satisfy 

H\ip) = 0, which for its coefficients translates into the difference equation 

\A f (jjl + AF(jj))\F(jjl + 4F( /U ))- 3 ^+4F M (0) - 2|A F (^)|F(^)- 3 Vv(0) 

+\A F (fi - 4F(ji))\F(n ~ 4F(/x))- 3 Vv-4fm(0) = -64 v / 67r 7 3 / 2 G^ mattcr ( / u)^(0) (5) 

where Aj?(/x) := |/i+F(/i)| 3 / 2 — F(/i)| 3//2 . (The Wheeler-DeWitt equation in a particular 
factor ordering is reproduced in the continuum or large volume limit where fi ^> F(fi) at 
large \x [32].) 

While such a difference equation of varying step-size F{p) is difficult to solve and 
analyze, it may be transformed to an equidistant equation up to factor ordering. For this, 
we use ipfj.+F(n) = V'A+i — ^'(f i )i J ' + ' ' * where omitted terms are higher derivatives of 
■0/5 := by p>{[i>) '■= J M dp '/ 'F(u). All derivative terms can be shown to be corrections 
of higher order in h, such that one may use equidistant shift operators in jl up to quantum 
corrections which one can absorb in the ordering choice [26J. For specific refinement models 
of power- law form F(n) oc fi x , we have an equidistant equation in the variable jl oc 
which would be proportional to volume /x 3//2 for x = —1/2. 

In the region around /i ~ 0, which is a strong quantum region, the equidistant difference 
equation may differ significantly from the original, non-equidistant one. Moreover, at those 
small volumes it is difficult to find a specific form of refinement function F(n) which could 
reliably capture refinements of small lattices. A power law f(p) oc p x with x < can 
clearly not be used there. It is thus fortunate that the general difference equation (J5j) can 
be used to conclude general singularity avoidance: In a solution scheme one can simply 
step over the values of the wave function at n — where the classical singularity would be 
[21 [5]. The arguments given in [2] for constant step-size are unchanged for an equation of 
the form (J5J) provided that F(0) ^ as it is required for a well-defined difference equation. 
Thus, quantum evolution of the wave function continues beyond the classical singularity, 
which as this general statement has been shown also for anisotropic [HJ H3] arid spherically 
symmetric models [4"5] . 
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2.2 Harmonic cosmology 



We are interested in quantum corrections as they arise when an initial semiclassical state, 
which solves a difference equation such as (jHJ) becomes more quantum. We can therefore 
make several approximations by ignoring quantum effects which are not expected to be 
strong or which would not change the qualitative effect of quantum back-reaction. An 
example for the first effect is the inverse (a -3 ), which differs from a -3 at small scales. 
As such scales are not reached in what we consider herejfl we can assume a 3 a~ 3 ~ 1. 
An example for the second effect is factor ordering choices, which will change the precise 
values of quantum correction terms but not their qualitative behavior. Thus, our analysis 
is not sensitive to the fact that the ordering of the Hamiltonian constraint in loop quantum 
cosmology is not fully fixed. Taking a square root in ([3]), we thus obtain a Schrodinger 
equation 



q I 3 

P4fll>(p, <f>) = ih—ip(p,(j)) = ±y ■j^7~^f(p)- 1 fnn(f(p)c)p\il>(j),<f>) =: ±Hip(p,(f)) (6) 

whose main trace of the loop quantization is the occurrence of the sine0 

There is a specific factor ordering for which the system becomes exactly solvable [8]. 
Since there are no general restrictions on the ordering, we choose this one for our analysis. 
Again, we emphasize that conclusions drawn are reliable only if they are insensitive to 
the factor ordering. For this, we use our canonical variables {f(p)c, V) and introduce 
J := V exp(if(p)c) which depends on c only in a form allowed by the loop quantization. 
(Using V and J as basic variables implies that the phase-space appears as a cone since the 
mapping from (p, c) to (V, J) is not one-to-one at V = 0. But outside V = 0, which is our 
main interest, we obtain a correct description.) For operators, ordering J as indicated in 
its definition, we have the commutation relations 

[V, J} = hJ , [V, J f ] = -hJ ] , [j, J ] ] = -2hv - h 2 (7) 

of an sl(2,R) algebra. Most importantly, the Hamiltonian is a linear combination of these 
basic operators, H = —ji(J— jt), which guarantees solvability, i.e. the decoupling of 



3 Inverse scale factor corrections are not dominant in isotropic models for large matter content (although 
they can be larger than sometimes assumed; see the appendix of [26). They do become important in 
inhomogeneous situations where local energies smaller than the total matter must be used. 

4 Thc harmless-looking absolute value on the right hand side would be difficult to compute explicitly 
for operators, which would make them non-local in specific p- or c-representations. Here, however, this 
will not present a problem at all: We simply need to require an initial state to have a decomposition in 
ff-eigenstates where only states with positive eigenvalues occur. Since H is preserved, this property will 
continue to hold for evolved states. Imposing the absolute value is thus not an issue of evolution, but of 
initial states which can be much more easily dealt with (for Ap^ <C \p$\ which we are going to use). One 
may also worry that superpositions are not allowed of the right kind: positive- and negative frequency 
(i.e. .ff-eigenstates with positive and negative eigenvalue, respectively) instead of left- and right-moving. 
We do not allow superpositions of the first kind as per the condition on initial states. We do, however, 
allow superpositions of left- and right-moving states (corresponding to the i-choice in ((6|)). Note, however, 
that we are going to compute primarily expectation values and fluctuations, which are of interest only for 
individual wave packets in a superposition, not for the total superposition. 
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quantum variables from expectation values. (We are ignoring factors of AttG/3 and 1 — x 
which would simply rescale in our solutions. They can be reinstated simply by multiplying 
H with 2a/47tG/3(1 — x), which we will do in our effective Friedmann equation.) 

We are still dealing with a difference equation of the type (jSJ), but in a specific ordering 
which will allow us to get access to interesting properties of its solutions much more easily 
via expectation values and fluctuations. The specific difference equation reads 

dip 1 
-p^ = ih— = Hijj = _(/2^j_ 2 -(fi + 2)^ A+2 ) 

which is a first-order version of 

Rather than solving for wave functions and then computing their expectation values 
and fluctuations of V and J, we compute those quantities directly. Expectation values 
satisfy equations of motion 

>-M._i (<J> + <J1)) , i^f^ ft -^> (! , 

which are not coupled to fluctuations or higher moments and can be solved easily. For 
arbitrary states, we have 

(V)(<f>) = ±{Ae-* + Be*)-±h (9) 

(J}(0) = ^(Ae-t-Be^+iH (10) 

with two constants of integration A and B and H = (H). 

We have to impose reality conditions to ensure that the variable c appearing in the 
complex J is real. Classically, this is equivalent to J J = V 2 , which also has to apply 
for operators: JJ^ = V 2 . (But note JU ^ V 2 .) Upon taking an expectation value, this 
condition relates expectation values and second order quantum fluctuations: 

\(J)\ 2 - ((V) + W = (AV) 2 - Cjj + \h 2 (11) 

where the V^-fluctuation and J- J-covariance 

(AV) 2 = (V 2 ) - (V) 2 and Cj j = ^(J.P + .Pj)-(J)(.P) (12) 

appear. If a state is semiclassical at least once, (AV) 2 — Cjj is initially of the order h 
and remains so because it is preserved in time [IB]. (See also the explicit solutions ffTB]) 
and ([21} below.) Thus, inserting the solutions © and (HO) yields AB = H 2 + 0(h). This 
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implies that A and B must have the same sign, and we can define A/B =: e 2S to write the 



is no singularity of vanishing volume but instead a bounce at <fi = 5. Large H, i.e. much 
matter in the universe, ensures that other quantum corrections not included here do not 
destroy the bounce. For instance, corrections in a~ 3 or from the factor ordering are not 
relevant in this case. For small H, however, there is no guarantee for a semiclassical bounce 
as (AV) 2 — Cjj could possibly compensate the positive H 2 and make AB negative such 
that (V)(<p) would be sinh-like. (There may be quantizations in which all states bounce 
irrespective of whether they are semiclassical even at only one time. But this would depend 
sensitively on factor ordering choices, which we do not consider a reliable property.) 

2.2.1 Free effective Friedmann equation 

As with any linear system, one can easily derive precise effective equations for the harmonic 
model of cosmology. They have the advantage of being more intuitive to interpret and thus 
give a more direct handle on geometrical aspects which would be difficult to decipher from 
a wave function. To derive effective equations [UJ HH], one first defines the quantum 
Hamiltonian as the expectation value Hq = (H) in a general state, parametrized by 
its expectation values and moments. The quantum Hamiltonian is thus a function of 
the expectation values and all infinitely many moments, determining their equations of 
motion via Poisson brackets as they follow from quantum commutators. However, all 
infinitely many equations of motion for the quantum variables are in general coupled to 
each other and to the expectation values, which usually makes these equations impossible 
to solve. One thus has to look for suitable approximations based on a truncation of the 
coupled system to a finite size of equations and variables. The truncated equations are the 
effective equations, which describe properties of the quantum system in regimes where the 
approximations used are valid. 

In a linear system such as our solvable model for a quantum space-time, no truncation is 
required as the variables automatically decouple into sets of finitely many equations. Thus, 
the quantum Hamiltonian provides precise effective equations which in our case follow 
from H e s = (H) = — \i(J — J), only depending on expectation values. (In the context 
of the effective Friedmann equation we drop brackets denoting expectation values.) For 
expectation values, the equations (jSJ) provided before are thus precise effective equations. 
It is sometimes useful to formulate them as evolution equations in proper time rather 
than internal time 0, as this can then readily be compared with the classical Friedmann 
equation. To derive this, we use H e g = — za/47tG/3(1 — x)(J — J) = p^, where we now 
have reinstated the numerical factors, and write the equation of motion 




contribution — \K). Thus, there 




+ <x. 



Here, we have eliminated J + J using the reality condition (TTTT) which implies 



- A ((J + Jf + (i(J - J)) 2 ) = J J = (V + h/2) 2 + (AV) 2 - G JJ + j 
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and thus 

J+J . L ( J-J 




2V + H V \i(2V + h) 



a 



(13) 



with a small a := ((AV) 2 — G JJ + ft 2 /4) / (V + h/2) 2 as a relative quantum variables. (This 
relation reflects the identity cos(/(p)c) 2 + sin(/(p)c) 2 = 1, corrected by factor ordering 
terms which arise from the quantized J. The contribution to V can safely be ignored 
for our purposes. Notice that the derivation of o shows that its form depends on factor 
ordering choices in the Hamiltonian and basic variables.) 

To reformulate this as a Friedmann equation, we use the expression for V in terms 
of a of (pQ) for a refinement function f(p) = fop x - Proper time then enters through^ 
4> = {4>, -^matter} = a~ 3 p<f>, which implies the corrected Friedmann equation 

dV / I Y(V\ AttGpI f AnG 




a J \Z{1 — x) J \ V I 6 a? \ 3 a? J 

for a free, massless scalar. As one can easily see, the quantum correction can simply 
be formulated as a term quadratic in the energy density of the free scalar [49J (see also 
[50j l5Tl |3T]). The bounce manifests itself since a = is possible if pf ree = \v\/ a& — Pent = 
3a^ 2 ^ Ax /8nG /q7 2 • For x = —1/2 and fo ~ £p, this is a Planckian energy density, but can 
be smallei^l for larger /o or x > —1/2. 

For a general interpretation, one should note, however, that higher powers of p are not 
the primary correction in loop quantum cosmology. The reason for the new term is the 
higher curvature corrections from the presence of the sine in the effective Hamiltonian, 
which in this model can be reformulated as a simple correction of the energy dependence. 
This proves the correctness of the effective Friedmann equation, as realized in [8] , but only 
for this specific matter content. While it is tempting to use ( 1T4"|) for all matter contents, 
just replacing p 2 /2a 6 with p, this would overlook corrections arising from quantum back- 
reaction which must be present in general. We will come back to these corrections in the 
next section, where we derive correct effective equations in the presence of self-interactions. 



5 As we are not including quantum corrections in the matter Hamiltonian, this relation is identical to 
the classical one. 

6 Moreover, the critical density would be a-dependent for x ^ —1/2. In this case, it is not just the 
matter density which plays a role for the bounce but also the underlying quantum gravitational state and 
its refinement. There is no a priori reason why only the matter density should play a role for the bounce. 
In fact, it is the spatial discreteness of the quantum representation which implies a repulsive force, and 
thus the quantum gravity state is the primary reason for the bounce. In the language of refinement models, 
including the coordinate volume Vo, we have fQa 2+4x — f(p) 2 p ~ Af~ 2 ^ 3 p — (a 3 /AA) 2 / 3 which is the value 
of an elementary flux, or roughly the area L 2 of an elementary lattice plaquette. In terms of L 2 , we have 
the critical density p cl i t — Z/^ttG^L 2 for any x, which clearly shows that the critical density is determined 
by the underlying lattice. For x = —1/2, L happens to be independent of a, and thus the critical density 
is a constant. For x ^ —1/2, on the other hand, the elementary flux depends on the total spatial size in a 
way determined by the lattice refinement, thus also making the critical density a-dependent. Requiring L 
to be constant and of the order £p implies that x = —1/2 and makes the critical density Planckian, which 
produces the parameter choices of [7]. However, this requirement is not necessary. 
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2.2.2 Quantum variables 



The solution procedure can be repeated for higher moments of a state, for which the same 
effective Hamiltonian H c s = — |z((J) — (™)) provides equations of motion. Hamiltonian 
equations of motion, for instance for volume fluctuations (AV) 2 = (V 2 ) — (V) 2 , are then de- 
rived by means of Poisson brackets which follow from expectation values of commutators: 
using {(V 2 ), (J)} = (\V 2 J])/ih = -i(VJ + JV) and {(V) 2 ,(J)} = 2(V)([V, J\)/ih = 
—2i(V)(J) we have {(AV") 2 , (J)} = —2%Cyj. Other Poisson brackets can be derived anal- 
ogously (see the App. [A] and [T2] for explicit expressions), which results in equations of 
motion 



— (AV) 2 = -C V j-C v j , — (AJ) 2 = -2C VJ , — [AJ) 2 = -2C 



±C VJ = -\ { Ajf-\c„-W? , ±C VJ =-\(AJ?-\ CjJ -(W? 

= -Cyj-Cyj. (15) 
They are solved by [^6] 

i i 

(16) 
(17) 
(18) 
(19) 
(20) 
(21) 

with independent integration constants c\,...,Cq. Reality conditions are already imple- 
mented for real cj since this implies C v j = Cyj an d3 (AJ) 2 = (AJ) 2 . Solutions can 
be chosen to saturate uncertainty relations in order to determine properties of dynamical 
coherent states. This analysis has been performed in [15J [TO] in order to shed light on 
the possible asymmetries of volume fluctuations before and after the bounce. For large 
H, which is required for a massive universe, state properties such as c 3 /c 4 and c 5 /c 6 are 
extremely sensitive to initial values. Thus, from current knowledge of the classicality of 
the universe one can practically derive nothing about the precise fluctuations of the state 
before the big bang. As we will see below, these same parameters are also crucial for the 



(AV) 2 (0) = 


i 

2 


[c 3 e~ 


-20 


+ c 4 e 2 ">) 


^(ci + c 2 ) 


(AJ) 2 (0) = 


1 
2 


[c 3 e~ 


-2(f> 


+ c 4 e 2 ^) 


+ ^(3c 2 - ci) - i{c^ - c 6 e 


(AJ) 2 (0) = 


1 

2 


[c 3 e~ 


-20 


+ c 4 e 2 ^) 


+ ^(3c 2 - ci) + iic^ - c 6 e 


c v M = 


1 
2 


[c 3 e 


-20 


- c 4 e 2 ^) 




c v M = 


1 

2 


[c 3 e~ 


-20 


- c 4 e 2 ^) 




CjM = 


1 
2 


[c 3 e~ 


-20 


+ c 4 e 2 ^) 


+ ^(3ci-c 2 ) 



7 Note that it is not the square root A J for which reality conditions are imposed because they are based 

On ((Jt)*) - (Jt)2 = (( J)2) _ ( J)2. 
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behavior near the would-be bounce at = of an interacting system. We emphasize the 
role of H: a large H keeps the bounce away from small volume and its strong quantum 
effects, but heightens the sensitivity to small changes in an initial quantum state. Small H 
would reduce the sensitivity, but also the bounce volume and thus introduce more quantum 
effects in the evolution. When discussing interactions in what follows we will see further 
implications of a large value of H. 

2.3 Interactions 

We are now ready to introduce an interaction term via a potential W(4>) of the scalar, 
which destroys the exact solvability and implies quantum back-reaction. There are several 
difficulties, which complicate a complete analysis. All this has been discussed in [12], and 
we present a brief review here: First, <fi in general will no longer serve as a global internal 
time, such that one can only study patches of solutions in which </> would be a monotonic 
function of a coordinate time. Especially since we have to address a question about long 
evolution times, several patches are in general required. Secondly, the 0-Hamiltonian p^ is 
now time dependent, as it is a function of <fi through the potential. While 



-P* = ±H = ± ^\J$- ^HW(0 (22) 

can still be used as a classical Hamiltonian, upon quantization solutions to the equation 
—p^ip = ihdip/d<j) = ±Hip will no longer provide exact solutions of the original constraint 
equation ([3]) which is a second order differential equation in <ft: for solutions ip of the 
linear equation, we have p^ip = ^p^Hip = H 2 ip ± [p^, H\ip ^ H 2 ip. Finally, there are now 
couplings between all the quantum variables, which requires a detailed analysis of effective 
equations. 

For a sufficiently small and flat potential, all these difficulties can be evaded: (i) While 
there is no global time, a small flat potential allows long monotonic changes of <fi which 
one can analyze in one patch. Different patches can be combined, but this requires the 
consideration of general states as initial semiclassical states in the first patch will change 
in each patch and possibly in the patching process, (ii) The commutator [p^, H] which 
determines the error one makes by formulating the initial second order equation as 
a first order Schrodinger equation ([6]) is proportional to hW and thus small for a flat 
potential. Moreover, since it arises from a commutator, one may view this extra term 
as simply a quantization ambiguity, or correct it by adding an extra contribution to the 
Hamiltonian H used for the linear equation, (iii) Coupling terms can be handled by suitable 
approximations in effective equations. This will be the main focus here. 

2.3.1 Effective equations 

To analyze implications for the bounce, we have to change the classical Hamiltonian (f22j) in 
two ways: we use canonical variables for general x, taking into account different refinement 
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schemes of loop quantum cosmology, and then introduce the holonomy variable J and its 
complex conjugate. This gives (again dropping, for now, numerical factors) 



1 J — J T/S/C 1 -^) 

H = V a ( J ~ 1)2 ~ vV(1 ~ x)w ^) = -^r - % -j^t w ^ + ■ ■ ■ (23) 

in an expansion by the potential. For this non-linear classical Hamiltonian, the quantum 
Hamiltonian depends on all quantum variables. It can be derived from a formal Taylor 
expansion of the expression H((V) + (V — (V)), (J) + (J - (J)), (Jt) + (Jt - (Jt))) in 
V — (V), J — (J) and Jt — (J'). From the result one can derive Hamiltonian equations of 
motion via Poisson brackets between expectation values and quantum variables, a formal 
procedure which can be seen as a shortcut to computing expectation values of commutators 

([•,#]>■ 

In our case, we obtain the quantum Hamiltonian 

7_J T/3/(l-x) 

3 2 4- r lAi+^X 1 -*') 3? U( 2 +a)/(l-aO 

-2» (1 -^ J - J + ~ °^ W ^ 

y3/(l-x) 

- 4 _((AJ) 2 - 2Cjj + (AJ) 2 )H/(0) + • • • 

where the omitted terms now also include higher moments. We will regard this truncation 
as our perturbative effective Hamiltonian, which is valid as long as higher moment terms 
are sub dominant. 

Due to the new coupling terms, equations of motion become lengthy and we do not 
present all of them here. Complete equations to the order used here and details of their 
derivation can be found in [12J for x = 0. The main equation of interest here is that 
for (V) because it is the one giving us information about the bounce. From the effective 
Hamiltonian, we find its equation of motion 

d(£> = _(J) + (jt) < J) + < Jt) { y ?l ^ )w(4>) (25) 

d4> 2 ((J)_(Jt))2 N ' W ^ I 

+:i ( i- W ( $ - £ W" +2 * )/(1 -" WffM 

+3 < ^ < ^ t> (l>) 3/ "-'((AJ)' + (A!)* - 2C J} )W(4,) 

~ 1 6 ^ + ^J V) {2+X)I(1 - X \CVJ ~ Cyj)W{<P) 
1 - X « J) - ( Jt})3 

O/T/^S/C 1 "^) 

: ((AJ) 2 -(AJ) 2 )H/(</>) 



3 /W^Va-^) 

+^ ^4 ; (Cyj + C V j)W{<t>) 

l-X ((J) - (Jt))2 V ™ VW 
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The evolution of the volume expectation value now depends on the behavior of the state 
via its moments. 

For semiclassical states, the correction terms are certainly small, but we recall that 
the singularity problem is a question about long-term evolution. Then, the state may 
change considerably and our equation, corresponding to the last phase before the big bang 
is reached, has to be applied to a state whose moments can be large. Additional correction 
terms will arise from higher moments neglected so far, but we can use the equation to test 
whether the fluctuations have an effect, and how this can occur precisely. Moreover, our 
conclusions can be extended to all orders in the potential and in quantum variables, which 
allows for arbitrary states [52]. We are here not interested in a precise analysis, which 
would seem premature, but in understanding the role of quantum corrections. Thus, we 
assume them to be small and check the self-consistency of this assumption. 



2.3.2 Correlations 

Let us thus assume that our state has for some time evolved by the exact solutions of the 
solvable model, corrected only slightly by the classical potential term and quantum back- 
reaction. We may thus use the solutions (EJ) and f TTOT) for expectation values and (|T6l) - (l2"TT) 
for fluctuations and covariances as zeroth order solutions to estimate the magnitude of 
quantum back-reaction. In particular, we look at where the free solutions would have their 
bounce, thus implying (J) + ( Jn) « 0. This already eliminates several of the terms in (1231) . 
in particular the classical interaction. However, the last two quantum correction terms 
remain, which further using (V) ~ H ~ ~\i{{J) — (J )) gives 



d0 ~ tf-a./a-*) V2 MA/(W + ^7^) Re ^(0) ) W&) . (26) 

This vanishes around ~ for background solutions satisfying C3 ~ C4, C5 ~ eg, which 
corresponds to unsqueezed states. For squeezed states, however, the time derivative of 
(V) does not vanish for the interacting system at the bounce of the free solutions. The 
interacting system cannot bounce where the free system would bounce. Thus, quantum 
back-reaction does affect the bounce for squeezed states. Whether or not and where the 
bounce may happen depends on quantum properties of the state, not just on expectation 
values as in the free case. Even for free solutions, there is not much control over the key 
parameters C3/C4 and c^/cq due to cosmic forgetfulness: it is precisely these parameters 
which determine the asymmetry of fluctuations discussed in [9j [10]. It will thus be even 
more complicated to constrain these parameters in an interacting system. To estimate 
the precise effect we have to know how squeezing develops for the interacting solutions; 
otherwise no reliable statement about the persistence of the bounce in the perturbed system 
can be made. 

In addition to cosmic forgetfulness, we have to consider the dynamical form of quantum 
variables and their coupling to expectation values in interacting states. Notice that the 
quantum corrections in (1231) split into two classes: the first three lines which include 
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( J) + (Jt) as a factor and which can be subsumed as changing the classical potential 
to W((f>)(l + 6%) with relative fluctuations 

e - o 2 + x ( A V ) 2 | o ( A J ) 2 ~ 2C ^~ + A J ) 2 6 ^ ~ 9x1 (27) 

and the last two lines which are present even when (J) + ( Jt) vanishes. 

These last terms are of special interest because they remain in the analysis of the bounce 
which led to 0261) . To see whether they are generically zero or, if not, how large they can 
grow, we need the equations of motion for Cyj and (AJ) 2 . This requires more Poisson 
brackets than presented so far, which are listed in App. |A]as recalled from [TJ]. In these 
relations, moments of third order appear (with triple superscript indices) which we ignore 
in our approximation together with the terms of order h 2 or higher. 

With these relations and the quantum Hamiltonian ( l24l) . we obtain 

-^(ReCVj) = ~(Re(AJ) 2 + Cjj + 2(AV) 2 ) (28) 
+2 <y )m ~ X) ( Cjj + ^^(AV) 2 ) W{cj>) 

+ ^(r^^ )(1+2 ' )/(1 ^' )(AV)2w/(0) 

m [vj_ w i mCVj ^(0) 

3 (VJ^V^^ + ^J)^ 



1-x ((J)-(Jt)) 2 
,(y)3/(i-*)((j) + (jt}) 



Im(AJ)W(0) 



and 



A(i m ( A J) 2 ) = -2ImCW + 4 ( 5 + X (f) 3/(1 7 } _ 3(2 + g) /^(i+toj/ft-,) | imC W{(j)) 
d0 l 1 KJ \l-x((J) - (Jt))2 2(l-x) 2W 1 



1 (I-*) 2 «J>-(Jt)W (J)-(Jt) 



-8*^4 s (2(AV) 2 + C Ji7 )W(0) 

-2z - 4—4^ ¥± Re(AJ) 2 W(4>) 

\l-x ((J)-(Jt)) 2 ; (J)-(Jt) 

6 <^^-- J (<i> + <i!>) Im(AJ) ^ ) . (29) 



i-a ((j)-(jt)) 2 
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Derivatives of the correlations by cannot vanish because they are determined by 
some of the fluctuations which, thanks to uncertainty relations, are non-zero. For instance, 
(AV) 2 appears in both equations and drives the change of squeezing. It is multiplied by 
the small potential, but also by factors of the volume. If we consider relative fluctuations 
(AV) 2 /V 2 of volume, which are nearly constant in each expanding or contracting branch 
of the free solutions, the dominating pre-factors are (V r )( 5 ~ 2z )/( 1_:r ') /((J) — ( J^)) 2 for ReCyj 
and (V r ) 3 ( 2 ~ a; )/( 1-:t: ) /((«/) — ( J" 1 ")) 3 for Im(AJ) 2 . Thus, especially at large volume, but also at 
the bounce itself, these terms can contribute to a significant change in time of squeezing. 
For x = —1/2, for instance, they behave as V A /H 2 and V 5 /H 3 , respectively. For the 
expanded Hamiltonian, on the other hand, only V 2 / H 2 W -C 1 is required, which allows 
the pre-factors to be large. Moreover, since the arguments can be extended to all orders in 
powers of the potential [52] , the condition on W to be small can be relaxed. When small 
volume is reached, squeezing parameters will be important in Eq. (1261) . 

Eq. (12"6"|) then shows that a large H has a protective effect for the bounce for x ^ 0, 
which is strongest for the limiting case x = —1/2. As before, we thus notice that a large H, 
which implies a bounce at large volume, preserves classicality better than small H, which 
is more pronounced the larger — x is. However, the influence of squeezing is always present 
and is required for a proper understanding of the bounce. Also this is weaker for x — — 1/2 
than for x = 0, but remains potentially large. We will now analyze the same effect from 
the point of view of an effective Friedmann equation which, unlike ( Tbll) . correctly captures 
state properties in the presence of interactions. 



2.3.3 Effective Friedmann equation 

Using the effective Hamiltonian, which is the truncation of the quantum Hamiltonian 
explicitly written out in ff2M . and the 0-evolution equation §Fo§ of (V), we can derive an 
effective Friedmann equation as before in Sec. 12.2. H . There are now several additional terms 
especially due to quantum variables, whose explicit dynamics we have not determined. We 
will therefore keep them in a general form as relative fluctuations such as fl2~Tj) and see how 
they affect the Friedmann equation. 

We first sketch the procedure and non-trivial changes due to the potential for the 
expanded Hamiltonian linear in W, but then use also the quadratic terms shown in App. [B] 
in the final effective Friedmann equation because they are expected if p 2 -corrections are 
present. In the effective Hamiltonian, it turns out that all relative quantum variables can 
be summarized by a single multiplicative correction to the classical potential: 

J - J 1/3/(1-*) 

with 

s = 3 2 + x (AV) 2 (AJ) 2 -2Cjj+(AJ) 2 _ 3 C VJ - C vj 

1 2(i-*) 2 {v y «i> - «it»2) i-x( V )({j)-{jt)) 

which differs from e\ in (J27J) by some coefficients. 
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To express dV/d<p in a similar form, we have to distinguish between two different 
contributions from quantum variables. One of them is E\, which, like 8\, just multiplies the 
potential. But there are additional contributions in the last two lines of (|25l) containing 
squeezing parameters, which are not of this form and which we write as^l 

(V)((J) - {J*)) i-* (V) 2 

We then have the evolution equation 

c\V 7 4- 7 7-1-7 „ ^ T/(4-a;)/(l-a;) 

% = -¥ + (jr^^-»^) ( i + fl) + v - TJ ^ rM . (32, 

With all numerical factors reinstated, we have 



H , s = 2 ^ (1 -J^l- ( ^Mi-iv y™ ww (1+g . 



v<$> 



which provides a quadratic equation for J — J in terms of p$. Its perturbative solution in 
W and fluctuations is 



and, via (U3l) . provides 



2^ = ± V 1 " ^^a^p + a + ({a°Wp?)>) 

with the energy density p = |a _6 p| + W(0) of the interacting scalar. 

The evolution equation f|5^|) now is to be squared for the effective Friedmann equation 
after replacing J + J. Here, terms quadratic in the potential become important, and so 
we should include them in the quantum Hamiltonian expanded to second order in the 
potential. The complete expressions are provided in App. [Bj which, with all numerical 



8 There is some arbitrariness in what one considers as relative fluctuations. Had we used (J) + (J^) 
instead of what appears in rji, then there would be a total factor of (J) + (J*) and r\\ could be subsumed 
in e%. It might, in fact, appear more natural to use, e.g., (Cyj + C v j)/(V)((J) + (J*)) instead of 
(Cy j + C v j) I (V) 2 . However, it is dangerous to do that because one would divide by a quantity, ( J) + (Jn), 
which then vanishes at the bounce. The quantity 771, on the other hand, simply absorbs factors which 
are already present in the equation or which do not vanish, and does so in a way which makes the final 
quantity dimensionless. Moreover, in each expanding or contracting branch (V) behaves approximately as 
(J) + (jt) but does not become zero. All relative quantum variables defined here are nearly constant for 
the free solutions in each branch, although they can change during the free bounce. 
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factors reinstated, provides the effective Friedmann equation 

p + e 1 W((f))-2^a t 



+ ^^fnl ■ (33) 

The ± corresponds to the sign of J+ J as it arises when the square root is taken. Coefficients 
denoted by Greek letters are relative quantum variables defined in the appendix or earlier 
in this section. From this equation, it does not follow in an obvious manner that there 
should be a bounce. As it stands, this equation may have no solution for a — or several 
ones. Also here, we clearly see the additional contribution by the squeezing term which is 
essential to decide about the bounce, i.e. to find out where, if at all, a — 0. If there would 
be no squeezing, i.e. 771 = = 772, then the scale factor would bounce for an energy density 
near the critical value of the free model, only slightly corrected by ei, 5\ and a. 

Instead of using the momentum and potential of 0, we can express the effective Fried- 
mann equation in terms of energy density and pressure P: 

a\ 2 8nG ( 1 , ^ Ap-P? 



a) 3 V 2 p + P 



x 771 



Pcrit V 2 4: p+P 

&(P-P) 2 ' 



2 p + P 
2itG (p - P) 



p 



Vt (34) 



where p crit = 3/87rG7 2 / 2 a 2+4x is the same as in the free case. However, p = p crit would 
give a bounce only in special cases. 
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3 Conclusions 



We have derived effective equations for isotropic models in loop quantum cosmology sourced 
by self-interacting or massive matter. In contrast to models studied so far in this context, 
this introduces quantum back-reaction which implies that the changing shape of a quantum 
state influences the motion of its expectation values. This effect, while initially negligible if 
one starts with a semiclassical state, is important for the long-term evolution of cosmology. 
The question of whether the classical singularity is replaced by a bounce can be reliably 
addressed only if quantum back-reaction is understood. Also the discussion of cyclic sce- 
narios for the universe, or any situation of long evolution times, requires these terms to be 
taken into account. 

We have derived perturbative effective equations which include quantum back-reaction 
terms, based on the availability of the exactly solvable model of [8] as free system. Com- 
pared to an effective Friedmann equation, which was previously available for models sourced 
by a free, massless scalar, there are new terms when the scalar has a potential. Phenomeno- 
logically, this implies that the effective Friedmann equation depends on the pressure of 
matter and not just on its energy density as in the classical or free quantum case. Com- 
pared to the free equation additional terms, which are independent of /o, do not arise from 
discreteness or repulsive gravity but from genuine quantum back-reaction. These terms 
do grow in long-term evolution and are especially important where the free system would 
bounce. We have proven that there is still a bounce provided that the state at small volume 
is uncorrelated, which presents a substantial generalization of what has been available so 
far. However, we have also demonstrated that the bounce is affected by quantum back- 
reaction of correlations, and much more work is required to see how generically a bounce is 
realized in this case. The now available equations, including those for moments, are suit- 
able for detailed numerical studies of the whole parameter space. However, the parameters 
relevant for the bounce are exactly the ones subject to cosmic forgetfulness [9j HQ] , such 
that a firm and robust statement about the bounce may be difficult to achieve. 

Moreover, our quantum back-reaction equations show that a generic state does develop 
strong quantum properties, especially correlations, which are important for the space-time 
picture at small volume and once a high-curvature region is traversed. This implies that 
space-time does not stay as classical as detailed investigations of the free model suggested; 
quantum geometry in general does become important. It makes itself noticeable not just in 
repulsive forces but also in genuine quantum variables such as fluctuations and correlations. 
While quantum evolution remains non-singular at a basic level of the whole wave function, 
in general there may not be a smooth geometrical picture such as a simple bounce of only 
the classical volume where genuine quantum variables would play no role. 

For isotropic models, there are several places where a large matter content H helps to 
reduce the effects of quantum physics, thus making the evolution appear more classical. 
(But note that the precise rate of reduction depends on the refinement scheme of loop 
quantum gravity which is only effectively included in isotropic models.) There are, how- 
ever, two important properties which prevent one from using this for a general physical 
conclusion to the extent that effects of quantum variables may be weak. First, already in 
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isotropic models, a large matter content enhances the sensitivity to quantum properties 
of an initial state [10] which, in long-term evolution, again makes quantum effects play a 
large role. Secondly, the scale of quantum effects is determined by the total matter content 
in the universe only for homogeneous models, while inhomogeneous situations are based 
on more local variables [17J. Those variables would generically be much smaller than the 
total matter content unless the spatial size of the universe is very tiny and comparable 
to the Planck volume. In either case, one cannot escape strong quantum effects when 
the small-volume behavior of cosmological models is considered. Thus, the behavior for 
small H is also of vital interest even though an isotropic model for a large universe would 
require large H. Similarly to a single massive partile in quantum mechanics, which does 
not show strong quantum properties in contrast to the microscopic dynamics of its ele- 
mentary constituents, quantum aspects of a large isotropic universe seem suppressed by 
its large matter content. However, an isotropic model does not fully capture the micro- 
scopic quantum gravitational dynamics, and it is this dynamics which must be shown to be 
non-singular. Within isotropic models one can obtain indications by looking at their prop- 
erties also for small H. While the basic removal of singularities based on an extendability 
of wave functions has been shown also for inhomogeneous systems of spherical symmetry 
|45j . which thus captures microscopic dynamics, this may in general not correspond to a 
smooth bounce. 

Acknowledgements: This work was supported in part by NSF grant PHY0653127. 

Appendices 

A Some useful Poisson brackets 

For equations of motion of quantum variables, we are using Poisson brackets 

{C VJ ,{V)} = iC V j , {C VJ ,(J)} = -i(AJ) 2 , {C VJ ,(.P)} = iCjj + 2i(AV) 2 

{C VJ ,(AV) 2 } = 2iG VVJ + 2i(J)(AV) 2 - \ih 2 (J) 

o 



{C V j, C V j} = AiG vvv + 6i((V) + \h)(AV) 2 - i(J)C vJ - i{P)C VJ + -th 2 (V ) + -th 3 



{C VJ ,(AJ) 2 } = -2iG JJJ -2i(J)(AJ) 2 

{C vj ,Cjj} = 2iG VVJ + 2i{(V) + \h)C V j-i(j)Cjj + i(P){AJ) 2 - l -ih 2 {J) 
{C V j,(AJ) 2 } = 6iG VVJ ' + %i((V) + \h)C VJ -2i(J)(AJ) 2 -ih 2 (P) 

for Cyj and 

{(AJ) 2 ,(V)} = 2i{AJ) 2 , {(AJ) 2 ,(J)} = , {{AJ) 2 ,(P)} = AiC VJ 
{(AJ) 2 , (AV) 2 } = MG VJJ + Ai(j)C VJ 
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{(AJ) 2 ,C V j} = 2iG JJJ + 2t(J){AJ) 2 

{(AJ) 2 ,C VJ -} = QiG VVJ + 8i{(y) + \h)C VJ -2i{P){AJ) 2 -ih 2 {j) 

{(AJ) 2 ,Cjj} = AiG VJJ + Ai{(V) + \h){AJ) 2 

{(AJ) 2 ,(AJ) 2 } = 8tG vvv + m((V} + lh)(AV) 2 + 8z((V} + ±h)Cjj-8z(j}C v j 

-8i{.P)c VJ + m 2 (v) + 2ih 3 

for (A J) 2 . They can be derived from expectation values of quantum commutators as 
discussed in more detail in [12]. 



B Expansion to second order in the potential 

To second order in the potential, we have the classical Hamiltonian 



H = \j~\{J- J) 2 -VW~*)W{(j)) 

J- J V 3 /^ 1/6/(l-x) 

= i ^-W + i ^-W 2 + --- 

2% (J- J) 2 (J - J) 3 

which implies the quantum Hamiltonian 

/ 1/6/(1-2) 5 _|_ r y2(2+x)/(l-x) 

"» = + w ' (fop + ^-z/(j-Jr {AVf (35) 

IS? VC^X 1- *) T/6/(l-a;) 

<C V j ~ C v j) + 6zf_- ((AJ) 2 - 2Cjj + (AJf 



1-x (J-J) 4 v VJ ' (J-Jf 

J — J 1/3/(1-0;) T/6/a-x) 

= — -T3X (1 + ^ W " + ! (737F (1 + * )W ' 2 (36) 

where "lin" denotes terms linear in W already presented in the main text. In the last line, 
we have defined 

5 + x (AIQ 2 (AJ) 2 -2C JJ+ (AJ) 2 18 CW - C vj f37) 
2 ' (l-^) 2 ((^)-(^)) 2 l-*(V>)«J)-</t» 

in addition to <5i in (1301) . 

From this Hamiltonian, we derive the equation of motion 

% = - £ T Z + (jr^v«/e-)(i + ei ) W -3^vW-»( 1 + ei)w . 

y(4-*)/(l-a!) T/(7-x)/(l-x) 

+ 1 7r7F" lH "- 3 1737F'' 2,1/ (38) 
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with 

5 + x (AVf | 10 ((AJ) 2 -2C JJ + AJ) 2 24 C VJ - C vj 

* := -4^-( AJ ) 2 + ^^ + ^ (40) 
(V)«J)-<Jt)) 1-x (V) 2 

in addition to ei in ( 1271) and 771 in ( 13TI) . 

The equality H e s = p^ now implies a fourth-order polynomial for J — J. To quadratic 
order in VF, its perturbative solution is 



.,. -r- . 1 + §-= — 5 — (1 + S 1 )W-^^—(l+25 1 + 5\-5 2 /2)W 2 J . (11) 



Inserting this into fl38|) and using (fl~3l) . we obtain 

U3/(l-i) T/6/ (1— a;) 
X I 1 + =— (1 + e 1 )W ^&W 2 



1 / p\ 1 y2(2+i)/(i-i) 

±o J 1 - 77? - v^+^/a-^Ci + - 5 + a 

2\ V 2 4 pi 



x 2 ^1 + 06 



2 



3- 



with 



1/6/(1-0=) 

+-^^nlw 2 (42) 



£1 := 1 + 25^ 52 + 51, (43) 

6 := ^1 + 6 1 - ^e 2 + Cl 5 2 - ^e 2 , (44) 
3 

£3 := ~Vi - + Vi^ - 2r] X 5i . (45) 



Reinstating numerical factors, this provides the effective Friedmann equation ( 1331) . 
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